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Abstract — The molecular analysis using the partial sequencing data of mt-DNA COI gene and nuclear 28S- 
rRNA gene is conducted to infer the phylogeney and geographical genetic variations in the genus Ryuthela. 
The genus Ryuthela is suggested to be the monophyletic and the most derived group in the subfamily 
Heptathelinae. Ryuthela iheyana, R. sasakii and R. ishigakiensis are monophyletic groups, respectively, but 
R. nishihirai seems to be paraphyletic to R. iheyana. Ryuthela spiders appear to be diverged into many phylo- 
groups with strong geographic associations. 
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Introduction 

Many members of araneomorph spiders are able to dis¬ 
perse great distances by aerial ballooning, but primitive 
mygalomorph spiders seldom disperse by ballooning 
(Greenstone et al. 1987). Therefore their interpopulation 
gene flow may be limited and prone to exhibit population 
divergence. DNA sequencing data showed that Aptosticus 
simus, the trap door spider with a limited dispersal capabil¬ 
ity, is geographically subdivided and divergent in the ab¬ 
sence of morphological differentiation (Bond et al. 2001). 
Furthermore, in the trapdoor spider species Antrodiaetus 
unicolor, species-level paraphyly was detected by DNA se¬ 
quencing data (Hendrixson & Bond 2005). 

Although I recognized four species in the genus Ryuthela 
by the male palpal morphology (Tanikawa 2013), species 
concept based on morphological distinctiveness appears to 
underestimate true diversity for spiders with limited disper¬ 
sal capabilities (Bond et. al. 2001). Because the genus 
Ryuthela is trap door spiders with a limited ability to dis¬ 
perse and have only a slight morphological differentiation 
among species (Tanikawa 2013), as in Aptosticus or 
Antrodiaetus, geographical divergence and/or species level 
paraphyly or polypohyly are expected. I conducted molecu¬ 
lar analysis using nuclear 28S-rRNA gene and mitochon¬ 
drial COI gene to infer the phylogenetic relationships 
among the species and to assess intraspecific and 
interspecific genetic variations in the genus Ryuthela. 

Materials and methods 

Sampling. The specimens used in this study were col¬ 
lected from throughout the distribution range of Ryuthela, 
that is, from Okinawa)ima Is. to Iriomotejima Is. of Ryukyu 


Isis., the southwest Japan (Fig. 1). Female and juvenile 
specimens were preserved in 99.5% ethanol at 4°C. Male 
adult spiders were preserved in 75% ethanol at room tem¬ 
perature, except for dissected right fourth leg preserved in 
99.5% ethanol at 4°C. Young males were brought back 
alive to the laboratory and reared until becoming adults. 
Two specimens of Liphistius sp. from Thailand and 
Malaysia, three specimens of Heptathela hangzhouensis 
from China, and two specimens of Heptathela yanbaruensis 
from Okinawa)ima Is. were used as the out group in 
phylogenetic analysis. Sampling data of the specimens used 
in this study are shown in Appendix. 

DNA extraction, polymerase chain reaction and se¬ 
quencing. Genomic DNA was extracted from muscle of 
legs of large individuals or the whole cephalothorax of 
small individuals using DNeasy Blood & Tissue kit 
(Qiagen, Inc., Germantown, MD). The mitochondrial 
cytochrome oxidase subunit I (mt-COI) partial sequences 
and nuclei 28S-rRNA (28S) partial sequences were used for 
phylogenetic analysis. Mt-COI was amplified using the 
primer combination CBl: 5' ^ TAT GTA CTA CCA TGA 
GGA CAA ATA — 3' (Jermiin & Crozier 1994) with 
HCOI-2198: 5' ^ TAA ACT TCA GGG TGA CCA AAA 
AAT CA — 3' (Folmer et al. 1994). The reactants were ini¬ 
tially denatured for 2 min at 90°C, proceeded with 40 cycles 
of 15 sec at 90°C, 20 sec at 50°C, 4 min at 72°C. 28S was 
amplified using the primer combination ZXl: 5^ — ACC 
CGC TGA ATT TAA GCA TAT ^ 3' with AS80P1: 5' 
^ AGA GCC AAT CCT TGT CCC GA ^ 3' (Bond & 
Hedin 2006). The reactants were initially denatured for 2 
min at 94°C, proceeded with 40 cycles of 30 sec at 94°C, 20 
sec at 50°C, 2 min at 72°C. PCR products were purified 
using the ExoSAP-IT (GE Healthcare Bio-Sciences, Co. 
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Fig. 1. Left: 50% majority rule eonsensus tree topologies obtained from Bayesian inferenee. Posterior probabilities are shown just after 
nodes. As for the details of topologies and distributions in eaeh island, see figs. 2-4 and appendix. Seale = 0.2 substitution/site. Right: 
Map showing the Ryukyn Isis. The underlined names show the islands from where Ryuthela speeimens were eolleeted. Seale = 100 km. 


Ltd., Buckinghamshire, England). The purified PCR prod¬ 
ucts were sequeneed using the BigDye terminator eyele se¬ 
quencing kit (ver.3.1) using the primer HCOI-2198 (mt- 
COI) or ZRl (28S) and analyzed on ABI 3100 or ABI 
3130x1 automated DNA sequencer (Applied Biosystems, 
Foster City, CA). Chromatograms were eheeked by eye 
using BioEdit Ver. 7.0.5.0.3 (Hall 1999) or MEGA version 
5.05 (Tamura et al. 2011). Sequenee alignments were done 
by MUSCLE (Edgar 2004) in MEGA. Overall mean p- 
distence was calcurated using MEGA by averaging the 
number of base differences per site over all sequence pairs. 

Phylogenetic analysis. The Perl script KAKUSAN 4 
(Tanabe 2011) and TREEFINDER (Jobb et al. 2004) were 
used to determine the appropriate model of DNA evolution 
by BIC for Bayesian analyses. MrBays ver. 3.1.2 (Ronquist 


& Huelsenbeek 2003) was employed to infer the phytogeny 
on combined data set. Four concurrent Markov Chain 
Monte Carlo (MCMC) ehains were run for 7,000,000 gen¬ 
erations, saving a tree every 100 generations. Topologies 
prior to In stabilization (“bum-in”) were discarded and pos¬ 
terior clade probabilities were computed from the remaining 
trees. The data of two Liphistius specimens were used for 
root estimation. 

Result 

In total, I sampled 270 speeimens from 75 loealities as 
shown in Appendix and obtained 578 bp of mt-COI, and 
697 bp of 28S-rRNA partial sequences from those speei¬ 
mens. About a half (9/20) of 28S sequences showed 
heterozygote. Overall mean p-distanees were 0.122 for mt- 
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Fig. 2. Left: Details of the topologies of elade A, R. nishihirai and R. iheyana. The symbols at nodes eorrespond to their loealities shown 
in right map. The OTUs are shown by the nemes of COI haplotype and 28S genotype. Seale = 0.1 substitution/site. Right: Map showing 
the loealities of speeimens belonging to eaeh sub-elade. Seale = 10 km. 


COI and 0.0154 for 28S-rRNA. The nucleotide sequence 
data are available in the DDBJ/EMBL/GenBank databases. 

The best-fit models of sequence evolution determined by 
KAKUSAN 4 were HKY85 + G for 28S, and for mt-COI, 
SYM + G for the first codon position, F81 + I for the second 
position and HKY85 + G for the third position, and gene 
proportional and codon proportional rather than other pat¬ 
terns of mixed model. Bayesian inference resulted in a 
phylogenetic tree shown in Figs. 1-4. 

Strongly supported (pp=100) three major clades were 
recognized in Ryuthela (Fig. 1), that is, elade A (comprising 
R. nishihirai and R. iheyana), elade B {R. sasakii), and elade 
C {R. ishigakiensis). Three species, R.iheyana, R. sasakii 
and R. ishigakiensis, were inferred to be monophyletic, but 
R. nishihirai seems to be paraphyletic to R. iheyana. The 


geographical population of R. nishihirai inferred to be a sis¬ 
ter of R. iheyana is living in the middle part of Okinawajima 
Is., far apart from Iheyajima Is. (Fig. 2). Ryuthela sasakii 
was not united with R. nishihirai inhabiting nearby but with 
R. ishigakiensis inhabiting in distant islands (Fig. 1). In 
each major elade, there were many minor clades arising 
mainly from mt-DNA variability with strong geographic as¬ 
sociations (Figs. 2-4), namely, Ryuthela populations were 
geographically divergent, though the relationships among 
them were ambiguous. 

Discussion 

The data set of the present study supports the monophyly 
of the genus Ryuthela, while the genus Heptathela is not 
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Fig. 3. Left: Detail of the topologies of elade B, R. sasakii. The symbols at nodes eorrespond to their loealities shown in right map. The 
OTUs are shown by the names of COI haplotype and 28S genotype. Seale = 0.1 substitution/site. Right: Map showing the loealities of 
speeimens belonging to eaeh sub-elade. Seales = 2 km, 10 km. 


monophyletic but is paraphyletic to Ryuthela. The 
phylogenetic position of Ryuthela inferred in this study is 
inconsistent with the cladogram inferred by Haupt using the 
morphological data (Haupt 2003, fig. 61). The present 
study shows that Ryuthela is a sister group of Japanese 
Heptathela and seems to be the most derived group in the 
subfamily Heptathelinae, while Haupfs (2003) inference 
showed that Chinese Heptathela (Sinothela in Haupt 2003) 
is a sister of Japanese Heptathela, and Ryuthela is a sister of 
all the other members of Heptathelinae (Heptathelidae in 
Haupt 2003). Haupt mentioned that ''Ryuthela possibly has 
hitherto unknown southern connection”, but the phylo¬ 
genetic position inferred in this study is not necessarily in 
such supposition but connection to the Eastern China along 
with Japanese Heptathela. Further examination should be 
done in future to clarify the phylogenetic structure of 
Heptathelinae and the position of Ryuthela, including the se¬ 
quencing data of Vietnamese Heptathela {Nanthela in 
Haupt 2003). 

There are many minor clades arising mainly from mt- 
DNA variability, which are distributed geographically 
nearby (Figs. 2-4), that is, Ryuthela spiders, at least 


females, are diverged into many geographical groups. This 
mt-DNA divergence supports the low mobility of female 
spiders as was observed in a previous study (Kikuya 1993). 
Although such a deep geographical divergence suggests a 
long history after expanding the distribution of Ryuthela, I 
could not separate them by morphology. There is a possi¬ 
bility of the existence of the cryptic species that are repro- 
ductively isolated with no clear morphological differences. 
However, as the males of Ryuthela wander to search for fe¬ 
males in the breeding season, nuclear DNA may exhibit a 
low divergence due to high mobility of males. It is there¬ 
fore necessary to analyze the population structure using nu¬ 
clear DNA markers, whose evolutionary rate are much 
higher than 28S-rRNA gene, to infer the actual amount of 
gene flow among geographical groups. 
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Appendix. Locality, haplotype, genotype and accession numbers of the specimens used for molecular work. 

Locality COI Accession 28S Accession No 

_ Latitude Longitude _ Haplotype _No_ Genotype _^ 

Ryuthela iheyana, Iheyajima Is. 


N27.04599 E127.95997 
N27.04599 E127.95997 
N27.04546 E127.95872 
N27.04178 E127.95813 
N27.04126 E127.96166 
N27.04126 E127.96166 
N27.04126 E127.96166 
N27.03361 E127.95542 
N27.03361 E127.95542 

Ryuthela nishihirai, Okinawajima Is. 
N26.73382 E128.22142 
N26.72760 E128.21060 
N26.68187 E127.89890 
N26.68187 E127.89890 
N26.68187 E127.89890 
N26.68187 E127.89890 
N26.67370 E128.18543 
N26.62611 E127.99059 
N26.62611 E127.99059 
N26.62611 E127.99059 
N26.62398 E128.06210 
N26.62398 E128.06210 
N26.61970 E127.91851 
N26.61970 E127.91851 
N26.60621 E128.14336 
N26.60621 E128.14336 
N26.60621 E128.14336 
N26.58921 E127.99792 
N26.58921 E127.99792 
N26.58921 E127.99792 
N26.58921 E127.99792 
N26.57086 E128.12697 
N26.57086 E128.12697 
N26.57086 E128.12697 
N26.56284 E127.98829 
N26.56284 E127.98829 
N26.55508 E128.05146 
N26.55508 E128.05146 
N26.54453 E128.08456 
N26.54453 E128.08456 
N26.52048 E128.01381 
N26.52048 E128.01381 
N26.52048 E128.01381 
N26.50046 E127.91703 
N26.50046 E127.91703 
N26.48550 E127.98256 
N26.48550 E127.98256 
N26.48550 E127.98256 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46752 E127.83551 
N26.46714 E127.84935 
N26.42298 E127.78798 
N26.42298 E127.78798 
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